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^ ■ ABSTRACT 

^ . We have made ROSAT VSVG and HRI X-ray observations to study the intraclus- 

' ter gas surrounding the powerful radio source Hercules A. The cluster is luminous in 

\ X-rays (Lboi = 4.8 x 10'^^ W), although apparently poor in optical galaxies, and the 

■ host of the radio source is the central dominating galaxy of the cluster. 

' The azimuthally-averaged X-ray surface brightness profile is well fitted by a mod- 

O ! ified King (/3) model, with core radius = f 2f ± 10 kpc and (3 = 0.74 ± 0.03, but 

■ the cluster is elongated parallel to the radio source, especially on the scale of the ra- 
dio lobes, and fits to individual quadrants give a core radius 50 per cent larger along 
the radio axis. Part of this elongation appears to be associated with enhanced X-ray 

^ ( emission superimposed on the outer radio lobes, which extend to just over 2rc- There 

are no obvious depressions in the X-ray emission coincident with the radio lobes, as 
expected if the relativistic plasma displaces the ICM. However, we show that these 
depressions may be quite weak, essentially because the main part of the lobes are 
^ , outside the cluster core. From the surface brightness profile for the PSPC data the 

X-ray emission extends out to ^ 2.2 Mpc radius. 

In the absence of the powerful jets (which must be a transient phenomenon on 
cosmological timescales), we would expect a cooling flow at the centre of the cluster; 
• but currently it must be substantially disturbed by the expansion of the radio lobes. 

\ The PSPC spectrum reveals a cool component of the ICM with 0.5<fcr<l keV in 

addition to the « 4 keV component detected by ASCA and BeppoSAX. The central 
cooling time could be as low as 2 Gyr if the cool component is centrally concentrated, 
otherwise it is around 6 Gyr. Cooling is significant on a Hubble time to a radius of 
about 90 kpc. The modelled central electron density of no = 1.0 x 10'' m^'^ is typical 
for modest cooling flows. 

Finally, we have detected faint X-ray emission from a compact central source, with 
size < 15 kpc and luminosity « 2 x lO'^^ W. 
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1 INTRODUCTION 

Hercules A (3C 348) is the fourth brightest DRAGN^ in the 
sky at low frequencies. At a low redshift of z = 0.154, its 
power at 178 MHz is Pitsmhz = 1.9 x 10^'^ W Hz"^ sr~\ 
with Ho ~ 65 kms^^ Mpc^^ and qo = (used throughout 
this paper). 



^ Double Radiosource Associated with Galactic Nu- 
cleus; sec Leahy 1 1993) or Leahy, Bridle & Strom 
I http://www.jb. man. ac.uk/atlas/ I for a full definition. 



Hercules A is identified with a very elong ated cD galaxy 
(e.g. ISadun fc Have^f 1993: B aum et alJll996l). with absolute 
magnitude —23.75 in the R-band JOwen fc La'in3ll989D . It 
appears to lie at the center of a poor, faint cluster, al- 
though measurements of the clu ster richness m ay be un- 
certain (Owen fc Laing 1989. : ..Yates. Miller fc Peacock.l98gl : 
iBarthel fc Arnaudlll99^ 

Radio galaxies with the radio and optical luminosity of 
Her A ne arly all show hotspot-d ominated Fanaroff-Riley II 
struc ture iLedlow fc Owedll99d) . but Her A IS IS an excep- 
tion (iDreher fc Feigelsonlll983r Its structure is dominated 
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by its twin jets, which are quite different in appearance. The 
eastern one is the brightest (highest flux density) radio jet 
in the sky, and contributes a substantial fraction of the ra- 
dio luminosity of Her A; a weaker jet to the west leads to 
a striking series of shells which dominate the western lobe. 
There are no compact hotspots. The structure is formally 
FR class I, but does not resemble typical FRI objects in 
detail. 

X-ray emission from the Hercules A cluster was first 
detected by the Einstein Obs ervatory iFeigelson fc Ber'3 
Il983l: iDreher fc Feigelsonlll984^ . They estimated the 0.2- 
4 keV luminosity as 2 x 1 0^^ W. This is typical of a rich- 
ness to 1 Abell cluster |^^ei|^5^. In their deprojec- 
tion analysis of the Einstein data on clusters of galaxies, 
IWhite. Jones fc FormanI lll997f) list Hercules A as a possible 
laxge cooling flow, although their best-fit model had a zero 
infl ow rate due to their low s patial resolution. 

iBarthel fc Arnaridl lll996l) have shown that Hercules A 
and other radio galaxies embedded in dense cluster gas are 
anomalously radio- loud (for their far-infrared luminosity), 
as expected because their confinement is more effective, re- 
ducing adiabatic expansion. 

In this first paper of a series of three, we report ob- 
servations of Hercules A i n X-rays made with th e ROSAT 
PSPC and HRI detectors. iGizani fc Leahvl ll200l hereafter 
Paper II) present new VLA observations of this powerful 
DRAGN, and we discuss the results derived from the spec- 
tral index and projected magnetic field/fractional polar- 
ization images all at l".4 resolution. In Paper HI (Gizani, 
Leahy & Garrington in preparation) we combine the results 
of Papers I and II to study the magnetic field of the clus- 
ter gas surrounding the DRAGN. Preliminary reports have 
already been published by Gizani & Leahy (1996, 1999). 
Our observations h ave also been discussed, in less detail, by 
ISiebert. Kawai ~ B rinkmami (1999). 

Sections 2 and 3 describe the observations and data 
reduction for the PSPC and HRI respectively. In Section 
4 we analyse the radial surface brightness profile to infer 
the density distribution in the Her A cluster. In section 5 
we compare the radio and X-ray structures to search for 
detailed correspondences, and review observations of related 
objects to put our results into context. Section 6 gives our 
conclusions. 



2 PSPC OBSERVATIONS AND DATA 
REDUCTION 

The total observing time with the PSPC detector was 8065 s, 
giving a live time of 7876 s. There were 5 observation inter- 
vals during 1993 August 19-31. These data were processed 
with the ASTERix package, which is the U.K. X-ray data 
analysis package running within the Starlink ADAM/ICL 
environment. 

The raw data, in MPE format, were pre-processed and 
converted to the ASTERIX internal format. Then t hey were 
'clean ed' using standard houskeeping parameters (lOsbornd 
Il993l) : to ensure good aspect solutions, we only accepted 
times with ASP_ERR< 1, which removed 3 seconds of data. 
We have found no significant difference by applying screen- 
ing to the spectral analysis for low particle background using 



the master veto rate (EE-MV) < 170. Therefore in the fol- 
lowing analysis we have included all data. 

2.1 Spectral Analysis 

The X-ray emission from Her A is sufficiently compact that 
the radial variation in sensitivity in the ROSAT field can be 
neglected. We therefore extracted a background-subtracted 
spectrum in ASTERIX using the normal 'point source' proce- 
dure. We used a source region with radius 4.7 arcmin, and 
a background taken from an annulus extending to 9 arcmin. 
The background-subtracted spectrum was corrected for ex- 
posure, dead time, vignetting, and for photons scattered out 
of the detect cell (assuming a point source model). 

The spectrum was processed with the f tools pro- 
gram GRPPHA using the binning recommended bv lBriel et alJ 
(1996). The data and response function were then read into 
xspec for spectral fitting. 

We fitted a MEKAL model jKaastral Il992^ 
with neu tral absorption using the xspec 11.0 cross- 
sections fBalucinska-Churc h fc McCammonI 119921) . and 
lAnders fc Grevesse (.1989) relative abundances. We do not 
consider Nh < 5.5 x 10'^° cm~^ as this is ruled out by 
Hi absorption llColgan. Salpeter fc TerzianI Il988h . which 
detects only cold Hi clouds and so is a l ower limit. We 
use Nh = 6.2 x 10^° cm"^ IStark et al.l (Il992l) which is 
in excellent agreement with the best fit Nh for our X-ray 
spectrum. A single-temperature fit for Z — 0.3, a metallicity 
value typical of clusters, which is consistent with the data, 
gives kT = 2.52l°i^ j^gY 

very similar to that found 
by ISiebert et al.l lll999l) . However, observations at higher 
energies give higher values for the best-fit temperatures: the 
LEGS detector on BeppoSAX gives ^ 3 keV in the 0.2-4 
keV band , while the 1.5-9.5 keV MEGS detector gives 
kT = 4.8 ± 0.6 keV and Z = 0.36 ± 0.15 llTrussoni et alJ 
[200 1). The 0.8-9 keV ASCA data of Siebert et al. are in 
excellent agreement with this, but Siebert et al. prefer a 
fit including a weak AGN component (consistent with the 
HRI image, see below), which gives kT = 4.25lJ;°^ keV 
and Z= 0.44 ± 0.13. At this metallicity the PSPG best-fit 
temperature becomes 2.7^q'4 keV, still not enough for 
reconciliation with ASCA & MEGS. Thus, as noted by 
Siebert et al., the combined results suggest that there is a 
multiphase cluster gas around Hercules A, and none of the 
single-temperature fits can be taken at face value. 

We have therefore tried a two-temperature fit to the 
PSPG spectrum with metallicity fixed at 0.44 and the hot- 
ter component fixed at 4.25 keV. For completeness we in- 
cluded a power-law component with F = 1.9 and normaliza- 
tion fixed to yield the point-source ffux derived in our image 
analysis (see below), although this is too weak (about 5 per 
cent of the total counts) to significantly affect the fit. With 
one extra free parameter, the fit improves by Ax^ ~ 3.5: 
the best fit temperature for the cool component is 0.71 keV, 
with a la range 0.44-1.01 keV; Fig.^shows this model over- 
laid on the observed spectrum for Hercules A. Although such 
a two-component model must be over-simplified, the PSPG 
response varies little with temperature or metallicity for gas 
above 3.5 keV, so the components plausibly stand for a ~4 
keV hot matrix, and a cooling-ffow region, which can clearly 
comprise a range of temperatures given the uncertainty in 
the best-fit value. The cool gas could also be associated with 
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Figure 1. The background-subtracted PI spectrum of Hercules A 
with the absorbed 2-temperature MEKAL model fitted. Model 
parameters are Nh = 6.2 X lO^" cm-^, kT = 4.25 and 0.69 keV, 
Z = 0.44 and z = 0.154, relative normalization 11.8:1 (corre- 
sponding to hard band iiux ratio 6.0:1). Reduced x^=l-21, with 
12 degrees of freedom. 



fc{T,Z,z) = 2.05 X 10"^" m^cts"Sr"\ 

The conversion factor between cts~^ and Wm~'^ (0.1-2.4 
keV, absorbed) is /b//c = 1.55 x 10^^* Wm-Vcts-\ and 
the cooling function, giving the conversion factor for 'bolo- 
metric' (0.01-100 keV) emissivity is: 

A(r, Z) = 1.90 X 10"^^ Wm^ 

Our model gives 

7/(ct s"^ sr~^) = 6.31 x 10~^£Af/(m"'' kpc), (1) 

where EM = J neiiH dl is the emission measure. For metal- 
licities in the observed range for clusters, 0.16 < Z < 0.5 
(e.g. Scharf & Mushotzky 1997), fc{T, Z) varies by less than 
10 per cent from our quoted value for 2.4 < kT < 10 keV. 
However f{T, Z) for the cool component increases rapidly 
with metallicity; for Z — 0.5 it reaches twice our quoted 
value for 0.5 < kT < 1 keV. If the cool gas is concentrated 
in the cluster core it likely has higher-than-average metal- 
licity and densities will be oner-estimated in that region by 
up to 50 per cent, but nearly correct at larger radii; if it is 
distributed throughout the cluster with a similar metallicity 
to the hot gas, the average density will be over-estimated by 
a few percent. 



the radio lobes, as seen in XM M and Chandra observation s 
of Vir A iBelsole et al."200lL Hyd A llNulsen et al.ll2002l) . 
and Per A ([Allen ct al. 2001a|), in each of which the radio 
jets seem to be dredging low-entropy gas out of the cluster 
centre. The excellent fits of single-temperature models to 
the MEGS and ASCA spectra do suggest that there is lit- 
tle gas at intermediate temperatures. The cool component 
contributes 11-21 per cent to the ROSAT count rate (ex- 
cluding the central compact source) with the higher contri- 
butions at the higher temperatures, and little dependence on 
metallicity. However, because it is dominated by line emis- 
sion, its normalization in terms of emission measure depends 
strongly on metallicity. 

The absorption-corrected bolometric luminosity of the 
Hercules A cluster is Lboi ^ 4.8 x lO^'' W (se e SectionlUl. The 
kT — L hoi relation for cooling flow clusters of lAllen fc Fabianl 
il99^ implies a temperature of around 3.6 keV for such a 
luminosity (with quite large uncertainty), consistent with 
the ROSAT, BeppoSAX and ASCA results. 

Unfortunately in our PSPC data there are too few 
counts outside the central resolution element to allow mean- 
ingful determination of the spectrum as a function of radius, 
so we cannot investigate any radial temperature gradients. 
Radial spectral analyses of the ASCA and BeppoSAX data 
by Siebert et al. and Trussoni et al. were also inconclusive, 
unsurprisingly in view of the lower resolution of these in- 
struments. We show below that the cool gas could represent 
the majority of the emission from the core of the cluster, 
but the phases could also be mixed, perhaps as a result of 
the expansion of the radio lobes. 

We derive the gas de nsity from radial profile fits to our 
X-ray images. Following iLeahv fc Gizanil ||200J) we evalu- 
ated the conversion factor between brightness and emission 
measure, for our assumed 'hot' component: for S in cts~^ 
in the PSPC hard band (PI channels 52-201): 



2.2 The PSPC Image 

Fig.Qshows that, due to Galactic neutral hydrogen absorp- 
tion, there are few source photons in the channels below 
0.5 keV; on the other hand the background is highest in 
these channels. Therefore the signal-to-noise ratio (SNR) is 
much higher in the hard band, and therefore we confine im- 
age analysis to this band. This has the added benefit that 
the PSF is smallest in this band and varies relatively little 
with energy. We created a background subtracted hard-band 
image using the same region as for the spectral analysis, and 
corrected for exposure, dead time, and vignetting. The po- 
sition scale was set by aligning the X-ray peak with the 
radio core. The offset of about 6 arcsec is consistent with 
the known pointing errors of ROSAT. 

The resulting image, smoothed with a gaussian 
of FWHM 20 arcsec, is presented in Fig. ^ Note 
that the fractional uncertainty in brightness in a 
photon image smoothed with a gau ssian is yj'l/2N 
(Hardcastlc, Worrall & Birkinshawl Il99^ . where A*' is the 
number of photons within the smoothing effective area 
(hereafter 'beam') A = 1.13(FWHM)^ In Fig. H A^ 0.9 
for the background, so the errors for faint emission are highly 
non- gaussian. 

2.3 The Radial Brightness Profile 

The faint outer halo of the Hercules A cluster is best seen in 
the radial brightness profile; however in this case the back- 
ground is detected with good signal-to-noise and it was nec- 
essary to take more care over background subtraction before 
we could obtain an accurate profile. In order to allow the 
background to be fitted as part of the overall profile fitting, 
our procedure was as follows. 

We located the clearly-detected point sources in the 40- 
arcmin diameter central portion of the PSPC field, using 
a simple peak finder applied to the image in Fig. H Peaks 
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Figure 2. The PSPC image of Hercules A. The image has been 
smoothed with a 20-arcsec FWHM gaussian. Resolution is 32 arc- 
sec FWHM, including the instrumental PSF. The image was pro- 
jected into B1950.0 coordinates in AlPS. The first contour is at 
5.29 X 10-1° W m ^ sr 1 , and contours are separated by ratios 

' sr^i = 0.9 ct 



of \/2. The background is 
beam" 1 . 



1.6 X 10^1° Wi 



and to the background model. Photons from an on-axis point 
source will not be affected by vignetting, so to compensate 
for the application of the correction to all the photons, we 
also applied a correction (normalized to unity at the centre) 
to our PSF model. 

There is a slope in the background across the PSPC 
field of view, with an amplitude of ~ ±15 per cent at 15 ar- 
cmin from the centre; however this is cancelled by azimuthal 
averaging. To get an idea of smaller scale fluctuations that 
might affect our profile, we summed the masked and cor- 
rected image between 0?15 and 0?3 in six radial bins; the 
rms between bins was 4 per cent of the (unsubtracted) back- 
ground, consistent with expected poisson fluctuations. We 
note that each radial bin has a similar solid angle to the sky 
region occupied by the cluster. Since it would be fortunate 
if the large-scale structure in the background were purely a 
linear slope, we conservatively assume there may be system- 
atic fluctuations in the azimuthally-averaged background at 
the 3 per cent level. 

The profile was binned to give binwidths > 5 arcsec, 
and a signal-to-noise ratio (SNR) of >5:1, subject to a max- 
imum bin width of 90 arcsec. The SNR was calculated after 
background subtraction, and the error used included our 3 
per cent systematic uncertainty, added in quadrature to the 
poisson uncertainties. All bins contained > 25 photons be- 
fore background subtraction. 

The binning and fitting procedure (Sectional was iter- 
ated a couple of times to ensure self-consistency. The final 
result is shown as the starred points in Fig. |4| 



within 5 arcmin of the cluster centre were not included in 
the list as they are likely to be substructure or random fluc- 
tuations in the cluster gas (the only clear point source in 
this region in our HRI image is not detected by the PSPC). 
We created a mask consisting of discs with radii ranging 
from 25 to 90 arcsec (depending on distance from the field 
centre, to allow for aberration) centred on each source. This 
was applied to the raw data, and a radial profile extending 
out to 0?3 was extracted, by integrating in rings around the 
cluster centroid using iradial, with a radial bin size of 2.5 
arcsec. 

The raw data must be corrected for vignetting, expo- 
sure and dead time. The correction as a function of radius 
is not explicitly available from ASTERIX, so we derived it by 
taking the ratio of profiles derived from the corrected im- 
age (from xrtcorr) and the image from xrtsub. To model 
the vignetting of the background, we also made a profile of 
the 'background' image created by xrtsub. (The same mask 
was applied in all cases). The model background includes a 
small (1.3 per cent) unvignetted contribution from the par- 
ticle background, but this negligible in our case. The model 
was derived from a 'spectral image' of the background an- 
nulus, in which the band was divided into 5 channels to ac- 
curately account for the energy-dependent radial vignetting. 
Since this region contains some source emission, we regarded 
the amplitude of the background as a free parameter in our 
fitting procedure (Section 0J. All profiles were converted to 
ascii files using ASt2txt, and further processed with our 
own software. 

The correction was applied to the raw data, to the as- 
sociated errors (initially assumed equal to poisson errors). 



3 HRI DATA ANALYSIS 

Our HRI observation was split into two periods separated by 
six months. The first, with a duration of only 1232 s, was on 
1996 February 18. The second, with the remaining 21516 s, 
contained 12 observational intervals during 1996 August 
27-30. Observations taken six months apart have opposing 
spacecraft orientations and frequently show large pointing 
offsets (D. Harris, priv. comm.). Since there were no strong 
sources in the field to allow accurate alignment, we decided 
to analyse only the second dataset. We first corrected it for 
the timing error in the standard analysis which affected the 
aspect solutions (using a script distributed by the Harvard- 
Smithsonian Centre for Astrophysics). The data was pro- 
cessed in ASTERIX and ftools. 

Beca use of the uncertain ty in spectral response of the 
HRI fe.g. lOreiner et al.ll999h . we did not attempt a spectral 
analysis. However, we did analyse the pulse-height analyser 
(PHA) distribution of source and background to try to op- 
timise the signal-to-noise ratio (SNR) . The background was 
determined in the annulus between 4 and 7 arcmin around 
Her A, which does not contain any background sources. The 
high HRI background makes it much less sensitive than the 
PSPC to very faint emission, so we decided to optimise the 
HRI image for the moderately bright emission around the 
radio lobes (at radii 50-100 arcsec), giving channels 3-7 as 
the best choice. These channels contain 83.3 per cent of the 
counts in a background-subtracted PHA spectrum of the 
source. 

We made an image with the full spatial resolution (0.5- 
arcsec pixels) using channels 3-7, subtracted the background 
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(0.00464 ct pixel"^, derived from the profile fitting in Sec- 
tional, s-'id converted to count rate by dividing by the expo- 
sure time. No vignetting corrections were applied but these 
are unimportant within 5 arcmin (the background is domi- 
nated by charged particles which suffer no vignetting). We 
also did not correct for detector quantum efficiency varia- 
tions, as these are not well determined at high resolution; 
we include a 4 per cent systematic error in the background 
to allow for this (c.f. Leahy & Gizani 200f ) . 

We converted from observed count rate to emission 
measure in the same way as for the PSPC data, except that 
we used the observed PHA spectrum to estimate the frac- 
tion of the IfRI counts that fall into our selected range of 
PHA channels. For our best fit model and our chosen PHA 
bands, the HRl is 2.86 times less sensitive than the PSPC, 
and this ratio varied by less than 2 per cent for any of the 
plausible models we tried. The lower sensitivity of the HRI 
is cancelled by the longer exposure time, so that the two 
images have nearly the same sensitivity to bright emission. 

Fig. 13 shows a greyscale of our A20-cm radio image 
(from Paper H) with contours from the HRl, with f 0-arcsec 
smoothing, superimposed. Positions were again set by align- 
ing with the radio; the HRI peak (after smoothing) was 
within 2 arcsec of the radio core. The HRI peak is not com- 
pletely symmetric, but we suspect that this is due to aspect 
errors during the observation. 

The radial profile for the HRI data was created using 
IRING in MPS as there were no spatially-varying corrections 
to be applied. The raw bin size was 1 arcsec and the profile 
extended to 250 arcsec. Full details of the procedure are de- 
scribed by Leahy & Gizani (2001). The profile was binned 
in the same way as for the PSPC, except that we assumed 
a 4 per cent systematic error in the background and used 
minimum and maximum bin widths of 2 and 50 arcsec re- 
spectively. The background-subtracted profile is plotted in 
Fig.|l] multiplied by 2.86 to put it on the same scale as the 
PSPC data. 

There are superficial differences between the PSPC im- 
age and the HRI image convolved to the same resolution, but 
these are approximately consistent with random noise: the 
apparent extension of the cluster peak emission, more or less 
perpendicular to the radio axis (see Fig. I^J is not confirmed 
by the HRI data. Whereas the PSPC seems to find a deficit 
of flux at r « 100 arcsec relative to a /3 model (see below), 
the HRI image shows a (marginally significant) plateau at 
about this radius, especially to the south-east, and in the 
'spur' extending to the north of the cluster, (which is con- 
siderably more diffuse in the smoothed HRI image than in 
Fig.HJ. 

To get a best-possible image of the cluster on arcminute 
scales, we made a weighted average of the PSPC and HRI 
images. The PSPC data were smoothed with a 20 arcsec 
gaussian as before, while the HRI image was smoothed with 
a 31.5 arcsec gaussian to give approximately the same res- 
olution, allowing for the « 25 arcsec FWHM intrinsic PSF 
of the PSPC. Both images were converted to instrument- 
independent brightness units. Weighting was by the inverse 
yj N/2 variance in the smoothed images, except that we did 
not let A'^ fall below the average background value. The 
larger smoothing of the HRI data makes the latter about 
twice as sensitive in the cluster core, while at larger radii 
the lower background of the PSPC means that the latter 
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Figure 4. The X-ray surface brightness profile of the Hercules 
A cluster from ROSAT PSPC (stars) and HRI (squares), after 
backgrounds have been subtracted. Error bars are la allowing 
for 4 per cent and 3 per cent uncertainty in the HRI and PSPC 
backgrounds respectively. 



dominates. Fig.|S]shows contours of this weighted image su- 
perposed on the same radio image as in Fig. |3 this time 
displayed with a logarithmic transfer function to emphasise 
the faint emission projected across the cluster core. 

The total radio emission extent is ~ 540 kpc, compa- 
rable to the size of the X-ray gas visible in the images, but 
the radial profile in Fig. 2] shows that faint X-ray emission 
extends out to at least 2.2 Mpc radius. 



4 MODEL FITTING OF THE SURFACE 
BRIGHTNESS PROFILE 

Allowing for their different resolutions, the HRI and PSPC 
profiles are in good agreement. The profiles differ within 
~ 20 arcsec of the centre due to the lower resolution of the 
PSPC, and are limited at large radii by the uncertainties in 
background level discussed in Sections 12. 2l fc |3 

We tried fitting the surface brightness profiles with var- 
ious models, to find the electron density distribution in 
the He rcules A cluster. Among these were a King or 13- 
model jKindll962l: ICavahere fc Fusco-Femianolll978l) . a de 
Vaucouleurs r^^'^ model and a combination (sum) of both. 
We confidently excluded the flat-topped 'King-law' profile, 
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Figure 3. A grey-scale image of tlie Hercules A radio emission at 1.5 GHz with contours of the HRI X-ray image superimposed. The 
radio image has resolution 1.4 arcsec FWHM, and the greyscale range is to 40 mjy beam~^. The HRI image has been smoothed with 
a gaussian of FWHM 10 arcsec, giving an effective resolution of 11.5 arcsec including the instrumental PSF. Contours are separated by 
factors of ^2 starting at 2.37 X 10"^ W m sr ^ = 3.07 ct beam ^ (3a above the background). The background level (which has been 
subtracted) is 1.63 X lO'^ Wm'^gj-i = 2.IO ct beam-^ 



as it failed to describe the brightness spike in the center 
of the cluster. The de Vaucouleurs model, as well as the 
model resulting from the sum of the latter with the King 
model, also failed, this time at larger radii. A somewhat 
better fit was obtained with an analytic 'cooling flow' model 
iBertschinger fc Meiksin|ll986h . but this still failed to rep- 
resent the central peak. We therefore concluded that this 
peak was a discrete feature and not just part of a contiim- 
ous brightness distribution. 

We obtained a good fit from a combination of a /3-model 
with a central point source. We emphasise that, although the 
/3-model was originally proposed as an approximation to the 
isothermal sphere, it provides an accurate em pirical fit to 
the e missivity in many 'cooling flow' clusters dWhite et al.l 
11997"). For the point source we used the theoretical HRI PSF 
i,David et al. 1998.1 . as none of the background sources are 
bright enough to give a measure of the PSF width. This does 
not give a perfect fit to the centre of the radial profile, not 
surprisingly, in view of the obvious elongation of the peak in 
the image. Following Leahy & Gizani (2001), we also tried 



fitting with a 'blurred' PSF to allow for aspect errors, but 
the improvement in the fit was only just significant at the 
5 per cent level, when fitting the HRI data alone; for the 
better-constrained joint fit to HRI and PSPC, the improve- 
ment was insignificant. Hence we use the nominal PSF in 
the following. The construction of the HRI PSF followed 
the method of Leahy & Gizani (2001), which allows for the 
finite pixel size of the original image. For the PSPC, where 
accuracy is less critical because the compact component is 
only weakly detected, we simply evaluated the PSF model of 
iHasinger et alJ il995l) at the centre of each 2.5-arcsec radial 
bin (assuming a photon energy of 1 keV). The results are 
given in Tabled Since the PSPC background model has a 
slight radial dependence, in Tabled we quote the value at 
the field centre. S is the point source flux and lo the central 
brightness of the /3-model. 

The PSPC profile fit was broadly consistent with that 
for the HRI, although the central peak is not clearly de- 
tected because the contrast against the cluster emission is 
much less with the lower resolution. To get a best estimate 
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Figure 5. A contour map of the summed ROSAT PSPC and HRI image of the Hercules A cluster in the 0.5-2 keV band, overlaid on a 
grey scale of the 20 cm radio map. The radio image is displayed with a logarithmic transfer function, between 0.3 and 300 mjy. Contours 
are separated by a ratio of y/2. The resolution of the X-ray contour map is 32", and of the radio image, 1".4. The first contour is at 
2.94 X 10"^'^ Wm~^ sr~^, and contours are separated by a ratio of \/2- The zero level has been subtracted. The boxed circle gives the 
smoothing beam of the HRI data, essentially the resolution of the X-ray map. 



Table 1. The final results from the model fit to the surface brightness profiles 





S 


h 




/3 




Background 




cts-l 


ct s"-*- sr"-*- 


kpc 






ct s"-*- sr"-*- 


HRI only 


0.0095 ± 0.0018 


9.6 ±0.9 X 10^ 


117 ± 19 


0.71 ± 0.08 


25.5 (25 df) 


1.05 X 10^ 


PSPC only 


0.0058 ± 0.0048 


10.7 ± 1.6 X 10^ 


112 ± 15 


0.72 ±0.04 


39.5 (27 df) 


9.59 X 10^ 


Combined fit 


0.0092 ± 0.0014 


9.6 ±0.6 X 105 


121 ± 10 


0.74 ± 0.03 


65.8 (56 df) 


9.60 X 10» 



HRI fluxes have been converted to the PSPC scale by multiplying by 2.86. 



of the cluster structure, we fitted both profiles simultane- 
ously, allowing the PSPC background to vary and keeping 
the HRI background fixed. The combined best fit is given 
in Table Q We use the parameters found from this in our 
further analysis. The best fit is shown plotted against the 
data from each instrument in Fig. |S| 



The combined fit is excellent for the HRI data. The 
HRI data contribute 25.6 to of the combined fit, with 
29 bins. The fit to the PSPC data is less good, = 40.2 
with 32 bins. Within a radius of 100 arcsec, the HRI data 
is as sensitive as the PSPC, and the good fit to the former 
suggests that the poor PSPC fit is just a random error. The 
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Figure 6. The best fit compact source + /3 model fit to the radial surface-brightness profiles of the Hercules A cluster. In each plot 
the combined model is given by the dashed line and the three components of the model (PSF, /3 model, and background) are given by 
dash-dot lines, (a) PSPC data, (b) HRI data. 



PSPC does reveal excess emission compared to the /3 model 
at 200-400 arcsec. 

The value found for P is within the range 0.5 
to 0.8 found for most cluste rs jjones fc Fo rmanI Il984t 
ICirimele. Nesci fc Treveslll997l) . ISiebert et alJ (Il999h find 
P = 0.63 and a core radius of 90 kpc (for our Ho), but this 
is based on fitting the HRI data only, and only in a sector 
away from the radio lobes. We replicate their result in Sec- 
tion |S3 The PSPC gives a much stronger constraint on /3, 
and at radii much larger than the radio lobes and likely un- 
affected by them, so it makes sense to average all azimuthal 
angles. 

To estimate the possible variability of the central 
source, we fitted the HRI and PSPC data separately, fixing 
the parameters of the j3 model and backgrounds to those 
found in the joint fit. The result is not well constrained be- 
cause the core is so weakly-detected by the PSPC; we find 
an insignificant increase of 13 ± 36 per cent between the 
PSPC and HRI observations. As expected from the small 
contribution of the core to the total flux, any power-law 
spectral component was undetected in the ASCA and Bep- 
poSAX spectra (Siebert et al., Trussoni et al.); as Trussoni 
et al. point out, this does confirm that the core is not heavily 
obscured. 

The integrated fiux in the PSPC profile within 500 arc- 
sec is 0.207 ± 0.005 ct s~^, which slightly exceeds the com- 
bined fiux in the model, due to the deviations from the model 
noted above. To get our best estimate for the total fiux of 
the cluster, we subtract the central source, apply a 3 per 
cent correction to infinite radius (found analytically from 
our P model fit), and an 11 per cent correction to the to- 
tal band (found by integrating the observed spectrum), PI 
channels 11-223, giving 0.227 ±0.006 ct s~\ corresponding 
to a 0.1-2.4 keV (absorbed) fiux density of 3.04 x lO"^'^ W 
m~^, a (rest frame, unabsorbed) luminosity of 3.00 x lO^'' W, 
and a bolometric luminosity of 4.84 x 10'^^ W (all calculated 
for our two-temperature spectral model). These figures are 



consistent with (but slightly lower than) the values derived 
from the much shorter RASS observation by Ebeling et al. 
(1998); they are in excellent agreement with the results from 
the long Emstem Observatory exposure analysed by White 
et al. (1997) (allowing for the different Ho,qo assumed in 
these studies). 

We saw above that relatively cool gas contributes 
/cool ~ 15 per cent of the cluster emission. This must have 
a density > 4 times higher than adjacent hot gas for pres- 
sure balance, and is also likely a more efficient emitter (i.e. 
higher f{T,Z,z)). Thus if it is radially segregated (as in 
cooling-fiow clusters), then it must nearly fill the cluster 
core; otherwise we would expect a two-scale core in emis- 
sion (the central source cannot be the cooling-fiow core, as 
its fiux is less than half that of the cool component). In fact 
11-21 per cent of the cluster emission comes from within 
0.65-0.89 core radii for P = 0.74. Thus complete radial seg- 
regation is barely consistent, even if /cooi is at the upper end 
of the range (which we recall implies (kT) closer to 1 than 
0.5 keV). 

We have found the electron density from the emission 
measure, EM, and the conversion factor derived in Sec- 
tion (recall this assumes kT — 4.25 keV). The electron 
density Ue of the P model cluster can be written as 



Taking — 1.2nH, we find a central electron density 
no = 1.0 X 10'* m^"^, suggesting a quite dense cluster, as typ- 
ical cluster central num ber densities are of order 10^ m~^ 
i Jones fc Formanlll984l) . If we have a cooling fiow, the cen- 
tral density is probably overestimated by 30-50 per cent, 
as noted earlier; but the densities derived from our model 
(with no as quoted) should be correct at r > rc, since the 
model correctly predicts the brightness at all radii, and the 
cool component makes at most a small contribution at large 
r. 
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Table 2. The results from the 3-component model fit to the on- 
lobe profiles 





rc 


/3 






kpc 






HRl only 


149 ± 19 


0.85 ±0.11 


28.5 (18 df) 


PSPC only 


148 ± 11 


0.84 ± 0.05 


23.3 (21 df) 



Table 3. The results from the 3-component model fit to the off- 
lobe profiles 













kpc 






HRl only 


98± 13 


0.65 ±0.06 


13.9 (13 df) 


PSPC only 


126 ± 11 


0.83 ±0.06 


39 (18 df) 



5 DISCUSSION 

5.1 Interaction of the radio lobes with the cluster 

gas 

We have looked for distortions in the cluster gas near the 
radio lobes, which might reveal interactions between the gas 
and the supposedly expanding lobes. As a first step we made 
radial profiles from the PSPC and the HRl data in four 
azimuthal sectors centered on the X-ray compact source. 
Two sectors were taken with position angles on the radio 
lobes ('on-lobe' profile), with a full width of 90° each, at PA 
98° and —82°. The other two sectors had position angles 
off the lobes ('off-lobe' profile), with the remaining position 
angles. Tables |5| and |3 show the results of the best compact 
source plus /3 model fit to the on-lobe and off-lobe profiles. 
The free parameters were 13 and the core radius Tc, the X- 
ray fiux from the compact source, the central brightness, and 
the background were kept fixed to the values found from the 
combined fit in Table Q 

Figs. Q and |S] show the on- and off- lobe profiles for the 
PSPC and HRl respectively. In each of the plots the com- 
bined model is given by dashed line and the three compo- 
nents (PSF, (3 model, and background) are given by dash-dot 
lines. 

Table H from the profil es on the radio lobes, shows that 
PSPC and HRl agree well, and show that the core of the gas 
distribution is significantly elongated in the direction of the 
lobes, as also visible in the images. There is less good agree- 
ment transverse to the lobe axis (Table suggesting that 
the profile departs significantly from a /3-model, especially at 
large radii where only the PSPC is sensitive. This is partly 
due to the faint excess emission to the north of the cluster 
visible in the PSPC image (Fig. |5J. Given the better fit of 
the model to the HRl data, we prefer the parameters derived 
from this fit for the inner regions of the cluster. If we take 
the off-lobe HRl fit as representing the undisturbed cluster, 
we can interpret the on-lobe fit as suggesting that the lobes 
have pushed some material out of the core, extending the 
apparent core radius and steepening the profile at slightly 
larger radii. 

We have seen that the X-ray emission is elongated along 
the radio axis. The maximum-resolution HRl image (Fig. 01 
shows discrete structures superimposed on the radio lobes. 



in particular arcs of emission perpendicular to the radio axis 
near the ends of each lobe. The brighter, western arc lies 
just beyond the brightest radio shell (or arc) embedded in 
the lobe. These features have been cited as evidence for a 
'radio/X-ray' interaction by Siebert et al. (1999). The ex- 
pected brightness from the smooth model of Table|21at radii 
74 and 88 arcsec (corresponding to the east and west X-ray 
arcs) is 2.2 ± 1.5 and 1.33 ± 1.3 HRl ct beam~^, whereas the 
observed peaks are 5.7 and 6.8 ct beam"'^. (Because Pois- 
son noise depends on the expected counts, errors are prop- 
erly attached to the model rather than to the data). At low 
count levels, positive deviations are more probable than for 
a Gaussian noise distribution, so the eastern 'arc' is not very 
significant, but the western arc seems to be real. The other 
structures noted by Siebert et al. are only marginally signif- 
icant, except for the most interesting, the ridge extending 
west from the core; unfortunately the latter is close enough 
to the core that it may be an artefact of pointing errors. 
As the discussion by Siebert et al. makes clear, there is no 
clear pattern in the relation between radio and X-ray fea- 
tures; e.g. the western jet is roughly superimposed on the 
X-ray ridge while the brighter eastern jet lies partly in an 
X-ray valley. Although Seibert et al. read the structure in 
terms of shells of material around the outer edges of the ra- 
dio lobes, the enhanced emission is almost all seen projected 
inside the outer boundary of the radio lobes, which makes 
interpretation problematic. If these features are analogous 
to the enhancemen ts around radio lob es seen in some other 
objec ts (e.g. Vir A. lBelsole et alj|200lt Hyd A. lNulsen et alJ 
[2002!), then we would expect them to be cooler than the 
ambient ICM. Indeed, they are plausible candidates for the 
spectrally-detected « 0.7 keV phase. 

Notably lacking in the X-ray images are any clear de- 
pressions in the cluster gas associated with the radio lobes, 
of the kind foun d in Cvg A JCarilh et aLllTgol I Wilson et alJ 
l2000l) . Per A iBohringer et alJ Il993l: lAllen et all 2001a^ 
and several other DRAGNs llMcNa,mara et al' l200d : 
iFinoeuenov fc JonedlioOll : lYoung et alJl2002l) . Given that 
the lobes of Her A are comparable in size to the cluster 
core, this is surprising, at least on the standard interpreta- 
tion that the lobes are relativistic 'bubbles' which displace 
the thermal gas. Certainly, we cannot interpret the radio 
bridges as components of a cylindrically symmetric bubble 
passing straight through the cluster centre; this puzzle is 
discussed in detail in Paper II. Here we focus on the 'bulbs', 
that is, on the bright, roughly elliptical outer lobes. Fig. |^ 
shows the predicted X-ray brightness from our /3-models for 
the Her A cluster, but with an ellipsoidal hole of zero X- 
ray emission matched in size to the radio bulbs. In Fig.|5Ji, 
the DRAGN is assumed to be in the plane of the sky. The 
brightness is reduced by a factor of almost two at the lobe 
centre. In Fig.|5j3, we assume the DRAGN is inclined at 50° 
to the line of sight, as estimated in Paper HI from our de- 
polarization data. In Fig. we take the cluster as elliptical, 
with the projected major and minor core radii as given by 
our sector fits. We assume the cluster is aligned with the ra- 
dio lobes in three dimensions, i.e. both have 50° inclination. 
Finally, in Fig. |5Jl the cluster elongation is in the plane of 
the sky, but the radio source remains at i = 50°. This gives 
very small distortions of the X-ray isophotes, although this 
geometry is contrived, implying as it does that the observed 
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Figure 7. PSPC data: Radial profiles in sectors centered on- (left panel) and ofli- (right panel) the radio lobes. 
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Figure 8. HRI data: Radial profiles in sectors centered on- (left panel) and off- (right panel) the radio lobes. 



alignment between cluster and DRAGN is merely a chance 
projection effect. 

The low signal-to-noise ratio in our HRI data, the low 
resolution of the PSPC, and the presence of the enhance- 
ments discussed above (which will partially mask any de- 
pressions), make it hard to exclude any of these models, 
although the case shown in Fig. is certainly disfavoured. 
Our HRI image (Fig. shows a ridge of emission along the 
northern edge of the western lobe, reminiscent of the generic 
features in Fig. |5| (but actually more prominent than in any 
of our models); but any ridges along the other lobe edges 
must be much weaker, at the level of the random fluctua- 
tions in the image. In short, our data are consistent with 
models in which the radio 'bulbs' are indeed empty of ther- 
mal gas, especially if we are correct in inferring that the 
radio axis is substantially out of the sky plane. Higher reso- 



lutions and deeper observations with Chandra should clarify 
the details of this interaction (Wilson et al., 2003, in prep.). 

5.2 Confinement of the radio lobes 

As the radio lobes of Hercules A lack anything resembling 
hotspots, it seems likely that they are expanding at around 
the sound speed, and so are confined by the thermal pressure 
of the ambient ICM, rather than by strong shocks. Therefore 
a measurement of the thermal pressure should be a mea- 
surement of the energy density of the lobes. The lobes are 
sharply bounded (e.g. Fig. 01, suggesting little entrainment, 
so the lobe energy should be dominated by relativistic par- 
ticles and fields. We therefore have the opportunity to test 
the minimum energy estimate Mmin with the more direct 
measurement from the external pressure. 
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Figure 9. Models of the X-ray emission from /3-model clusters with an ellipsoidal hole of zero emissivity, representing the radio lobe. In 
all cases the hole has major and minor projected axes 74 and 55 arcsec, its end is projected 96 arcsec from the cluster centre: the dashed 
ellipse shows its outline. The greyscale is linearly proportional to the model brightness, and the solid lines are X-ray isophotes separated 
by ratios of \/2. We take /3 = 0.74 for all the cluster models, (a) Spherical cluster with rc = 43 arcsec; radio axis inclination i = 90°. 
(b) as in (a), but i = 50°. (c) Elliptical cluster model with projected core radii 53 and 35 arcsec. Cluster and radio axis inclination both 
50° . (d) as in (c) but cluster assumed to be in the plane of the sky, while the radio axis remains at j = 50° . 



As in Section 15.11 we focus on the bulbs rather than 
the faint inner lobes. Although the ends of the bulbs are at 
almost exactly the same distance from the centre, the more 
elliptical western bulb reaches closer to the middle and so on 
average is closer to the centre. Following the methodology 
of Leahy & Gizani (2001) we have calculated the minimum 
pressure (Pmin = Umin/S) and the thermal pressure at the 
deprojected distance of the eastern and western bulbs of 
Hercules A, assuming i « 50°. For the calculation of the 
thermal pressure, Pth = l-Qn^kT, we have used the results 
of the combined fit of the surface brightness profile model 
(see Table The results are listed in Table 2] P is the 
deprojected distance of the bulb centres, while Bme is the 



Table 4. Minimum energy parameters and the thermal pressure 



Component 


R 


Bme 


p . 

mm 


Pth 


Bms 


Pms 




arcsec 


nT 


pPa 


pPa 


nT 


pPa 


East bulb 


84 


0.79 


0.19 


2.3 


1.99 


1.05 


West bulb 


71 


0.73 


0.17 


3.0 


1.86 


0.92 



minimum energy magnetic field. In addition to the conven- 
tional minimum energy parameters we also list the values 
Bms and Pms found using the methodology presented in 
Myers & Spangler (1985), that is, assuming the electron en- 
ergy spectrum extends down to a Lorentz factor of 10. 
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As we can see from Table |11 the radio lobes of Her- 
cules A would be significantly underpressured with the stan- 
dard minimum energy assumptions for the radio plasma: 
electron/positron composition (i.e. no protons), filling fac- 
tor unity and a low-energy cutoff corresponding to 10 MHz 
emission. Thus minimum energy is a severe underestimate 
of the actual energy content of the lobes, by an order of 
magnitude. 

Such a deficit in the minimum energy has been 
established for a long tim e for FRI DRAGNs (e.g. 
iFeretti. Perola fc Fantil[l992^ . but FRI jets and tails were 
often modelled as turbulent flows with substantial en- 
trainment, so the result could be attributed to mixture 
with the ambient gas and/or to a low filling factor for 
the radio plasma. This was ruled out in some FRIs by 
the the discovery of X-ray cavities fBohringer et al.l ll993l: 
LMcNaniara ct al. 2000), and by deficits in cases li ke Her A, 
where the lobes are sharply bo unded (iHardcastle fc WorralJ 
l2000al: iLeahv fc Gizanil I2OO2I) . This points to an error in 
the microphysical assu mptions, as discussed in detail by 
iLeahv &: Gizanil i200ll) for 3C 388. In that case we argued 
that the extra energy is mostly in the form of particles, and 
as a result the magnetic field should be below equiparti- 
tion and therefore does not play an important role in the 
lobe dynamics. The results of the Myers & Spangler calcu- 
lations show that an alternative option is available for Her A: 
equipartition can be (almost) preserved if the particle energy 
spectrum extends down to trans-relativistic energies. The 
difference arises essentially because the steeper spectrum of 
Her A gives a larger increase in energy on extrapolating to 
low Lorentz factors. 



5.3 Cooling and heating of the cluster gas 

The intracluster gas is densest in the core of the clus- 
ter and therefore the radiative cooling time icooi, due 
to the emission of the observed X-rays is shor test there 
jFabian. Nulsen fc Canizarejll99lt lFabian|[l99i) . For iso- 
barically cooling gas 



kT{nH +n^ + nHe) 5.75 kT 



A(r, Z)nenH 

-1 



A(r, Z) ne 



(3) 



10 



11 



103 m- 



T 



10«K 



yr 



where A(r, Z) is the cooling function (see Leahy & Gizani 
2001). The last form is the approximation of^razin ( 1986), 
which begins to become inaccurate for temperatures below 
3 keV. For the centre of the Hercules A cluster the cooling 
time is 6.4 Gyr if the hot phase dominates, but ~ 2 Gyr if 
the core is dominated by the cool phase. Since this is less 
than the age of the unive rse, this would normally imply the 
presence of a cooling flow ^ISarazinll986^ . The cooling radius, 
i.e. the radius at which the cooling time equals 10^" yr, is 
~ 90 kpc, somewhat less than the core radius, consistent 
with the relatively good fit to the /3 model. 

■Pore s ct al. (1998) find that the core radii of cooling 
flows clusters is generally smaller than in non-cooling fiow 
clusters, with the division around 130 kpc, which also sug- 
gests that Her A has a cooling flow. However, the Peres et 
al. core radii are fltted values for the (assumed isothermal) 



cluster mass distribution, and are not directly comparable 
to gas core radii derived from /3-model fits. 

It now appears that in many or all clusters, some agency 
prevents fully formed cooling flows from developing, so that 
little gas cools t o much less than half the ambient cluster 
temp erature fsee lMolendi fc Pizzolatoll200ll : |Peterson et alJ 
[2Q03 for recent analyses). Energy input fr om DRAGNs is 
a popular candidate for this agency, (e.g. *Churaz ov et alJ 
^002, and references therein). Attention has mostly been 
focussed on weaker, quasi-continuous AGN activity rather 
than major episodic outbursts as in Her A, but the la tter 
may also play an important role (e.g. lDavid et al]l200lf) . 

Although our results show that the normal state of the 
Hercules A cluster is a cooling flow, a region a few times the 
size of the core is now affected by the expansion of the radio 
lobes, and it is easy to see that the 'waste' energy supplied to 
the environment by the DRAGN currently far outweighs the 
energy loss from X-ray radiation. If we assume quasi-isobaric 
(and so subsonic) expansion, so that the cluster pressure 
distribution is not disturbed, the average rate of work done 
by the DRAGN is dW/dt = J^^^^^pdV/r where r is the 
lifetime. As in Section [5. II we focus on the outer 'bulbs' of 
the lobes, and we model these as prolate ellipsoids inclined 
at i ~ 50°; numerical integration gives a 'pressure content' 
of ~ 5 X lO'^^ J. Writing r = rg x 10** yr, this gives dW/dt ~ 
1.6rg~^ X 10^* W. The combined kinetic power of the jets K = 
AdW/dt, as it equals the enthalpy flux into the lobes (we 
assume that the lobes are dominated by relativistic plasma, 
so It = 3P). We can also parameterize the problem in terms 
of the radiative efficiency of the DRAGN. With the radio 
luminosity Lradio ~ 3.8 x 10^^ W, 

e = -/^radio/if ~ O.Oerg, 

and we can write dW/dt — I/radio/4e. Estimates for the 
lifetimes of powerful DRAGNs are typically in the range 
0.1 5;r8 ^ 1, corresponding to an efficiency of 0.6-6 per cent 
for Her A. Given that the radio luminosity of the DRAGN 
and the X-ray luminosity of the entire cluster are almost 
the same, the cluster gas is currently gaining energy at 3-30 
times its radiative loss rate. Revnolds. Heinz fc B egelmanI 
(j^OOl) point out that initially the work done may appear as 
potential energy (some of the gas is lifted out of the clus- 
ter core); in the long term much of this will be re-cycled to 
heat via turbulent dissipation. In supersonic expansion, of 
course, there will be at least some direct heating by shocks. 
In either case, any cooling of the cluster will be temporarily 
reversed. When the radio outburst is over, the cluster will 
return to its cooling state. At present, most of the mechan- 
ical energy is being deposited well outside the cluster core, 
and it is possible that the core is already cooling again, even 
while the outer part of the cluster is heated. 

If the efficiency is around 1 per cent, and the AGN is 
periodically active with a duty cycle of around 4 per cent, 
then on average radiative cooling of the cluster is balanced 
by mechanical heating. The cooling time of the cluster would 
then provide a rough measure of the time since the last out- 
burst. Such a scenario appears to suffer from fine tuniiig 
(Fabian ct al. 2001.) , but, as suggested by .Churazov ct alj 
Il2002f) . feedback via accretion onto the AGN from the cool- 
ing gas could naturally regulate the cycle of heating and 
cooling. 

Although it is most likely that Her A is expanding tran- 
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sonically, it is interesting to consider how compression of the 
cluster gas behind the bow shock of a DRAGN will affect the 
cooling time. We have found how the Mach number changes 
the cooling time, by using the Sarazin ( 1986) approximation, 
and the change in tem perature and density as a function of 
shock Mach number IILandaxi fc Lifshitzlll959ll . For shock 
normal Mach numbers M < 6.5 the cooling time is initially 
somewhat reduced, but only by up to around 40 per cent, 
which is not nearly enough to make radiative losses impor- 
tant on the timescale of the AGN ( < 10* yr). Over the longer 
term, the material will re-expand adiabatically to a condi- 
tion with a longer cooling time than before. Higher Mach 
numbers will increase the cooling time from the start, but 
for > 10 the shocked gas will be too hot to be visible in 
X-rays. We note that in the case of Hercules A the pre-shock 
central cooling time would have been < lO^'' yr whatever the 
Mach number, i.e. the apparent cooling flow cannot be an 
artefact of the DRAGN. 



5.4 The magnetic field in the cluster core 

In Paper HI we estimate the strength of the intracluster 
magnetic field as a function of radius, using our X-ray data 
and the Faraday dispersion from our radio data. In this sec- 
tion we attempt a cruder estimate, using the dispersion of 
the rotation measure ct rm in the central = 43.4 a rcsec. 
We follow the analysis bv lOarrington fc Gonwavl I^1991^ . who 
show that the expected aRM for a source behind a Faraday 
screen with a tangled magnetic field with rms B nT (inde- 
pendent of radius), is given by 



O'B.M 



0.0081 



(ne_B||) Idz 



1/2 



(4) 



where z is along the line of sight, I is the coherence length, 
(both in kpc), rie is electron density in m~^, and By is the 
component of the field parallel to z, equal to B/y/Z for 
isotropic tangling. For a source at the center of a cluster 
who se radial profi le is described by a /3-model this becomes 
fe.g. lFelten|[l996h 



0.0081B||nor'-''Z^/^ / 7rV^r(3/3 - 05) 
(l + rVr§)(6/5-i)/4 V 2T{W) 



(5) 



where no is the central electron density and T is the Gamma 
function. Figure IIUI is a histogram of the rotation measure 
for all significant pixels in Her A within the core radius, 
from our radio data (Paper III). We find (j_rm = 214 rad 
m~^, and the coherence length I near the centre is about 
5 arcsec or 14 kpc (Paper HI). We take a fiducial value 
of r = 80 kpc (28.6 arcsec), since most of the pixels with 
valid rotation measures are towards the outer part of the 
region. Then, using no, /3 and Vc from our combined HRI and 
PSPC fit, the estimated B ~ 0.2 nT. Hence the magnetic 
field in the central region of the cluster of Hercules A is 
lower t han in extreme coo ling flow clusters (e.g. the 3C 295 
cluster, lAUen et al.l|2001b^■ but within th e typical range for 
clusters in general dCariUi fc Tavloill2002fl . 



5.5 The central source 

As we have seen in Section 2] the improvement of fltting a 
blurred PSF to the central peak is insigniflcant. This means 
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Figure 10. Histogram of the rotation measure for all significant 
pixels within 43.4 arcsec of the radio core of Her A, from our radio 
data (Paper HI). 



that both a nominal PSF and a blurred PSF are both good 
fits of the weak central peak in the HRI image (for example 
see Fig. 1^1. AH we can say is that we have a central source 
with a count rate of 0.009 ct s~ in our hard band {Tahle^. 
The peak could be either an AGN or a galaxy-scale cooling 
fiow. We will discuss each case separately: 



5.5.1 The central source as an AGN 

We have already noted that there is no sign of an 
obscured AGN in t he high-energy BeppoSAX spectrum 
iTrussoni et al.l200ll) , so we assume a power-law model with 
only Galactic absorption. Assuming a typical photon index 
of r = 1.9, the unabsorbed 0.1-2.4 keV luminosity of the 
central source is then 2.0 X 10="^ W. At 5 GHz, the radio 
core luminosity is 7.1 x 10^^ W Hz~^. The point source ap- 
pears a few times over-luminous in X-rays relative to the 
correlation between Lx and radio core power for galaxies 
iSiebert et al.lll99d . corrected to our Hq), although it lies 
within the general scatter. We note that Siebert et al's cor- 
relation is formed mostly by broad-line objects, while Her 
A is a low excitation rad io galaxy (Tadhunter et al. 199^; 
Ijackson fc Rawlingj|l997^ , with a steep spectrum (optically 
thin) radio core at GHz frequencies (Paper II). Although X- 
ray variability is expected in an AGN we have not detected 
any during the interval of the two and a half years between 
our PSPC and HRI observations. AG N X-ray emission could 
be sync hrotron or inverse Compton. iHardcastle fc WorralJ 
i2000lJ^ have argued that inverse Compton is required, at 
least in FR Is, if the unified model of low-power radio galax- 
ies and BL Lacs is correct. 
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5.5.2 The central source as a cooling flow 

The extra broadening by the PSF that fits the data best 
corresponds to the nominal PSF convolved with a gaussian 
of FWHM 5.0 arcsec. Therefore the width of the cooling 
flow peak, if any, would have to be around 5.0 arcsec or 
less. In contrast the cluster core radius found from the fits is 
43 arcsec, an order of magnitude larger. This would give the 
cluster emission a 'prussian hat' look: essentially the central 
peak is on the scale of the central cD galaxy rather than the 
cluster. Small central peaks like this are sometimes seen in 
cooling flows. The X-ray luminosity is just greater than the 
largest value for an individual galaxy in the Einstein Ob- 
servatory sample of lEskridee. Fabbiano fc KimI 1^9^, but 
this sample excluded cluster-centre galaxies. This suggests 
that if the central source is thermal emission it relates to 
the cluster and/or the AGN, rather than solely to the host 
galaxy. In this case the central source, at around 5 per cent 
of the cluster luminosity, would presumably contribute to 
the cool component detected in the spectrum, but could not 
produce all of it, as previously noted. 

6 CONCLUSIONS 

Our ROSAT PSPC and HRI X-ray observations of the intr- 
acluster gas in the Hercules A cluster have revealed an ex- 
tended X-ray emission, extending to radius 2.2 Mpc (PSPC 
data) and a weak central peak (HRI data). 

Based on discrepancies between single-temperature fits 
to the PSPC, ASCA and BeppoSAX spectra, we argue that 
the intracluster gas contains at least two phases, domi- 
nated by a « 4 keV component (seen by ASCA and Bep- 
poSAX), but with around 15 per cent of the emission in the 
ROSAT band from material at 0.5-1 keV. This model sig- 
nificantly improves the PSPC fit compared to one with a 
single temperature. The 0.1-2.4 keV (absorbed) flux den- 
sity is 3.0 X 10~^^ W m~^, and the bolometric luminosity is 
4.8 X lO^'^ W. There is no evidence for absorption above the 
Galactic value. 

We have detected X-ray emission coming from a com- 
pact source at the centre of the cluster, with a 0.1-2.4 keV 
luminosity of 2 x 10^® W. Our data does not distinguish be- 
tween an AGN and galaxy-scale thermal emission (perhaps 
a central cooling-flow spike). The central electron density 
(excluding the central source) is no ~ 1.0 X 10" m-^ The 
central cooling time is in the range 2-6 Gyr, depending on 
whether or not the cool gas phase is concentrated in the 
centre. The cooling radius is 90 kpc. The cluster core should 
therefore be a cooling flow if we ignore the effect of the 
DRAGN. Indeed, most of the time, the DRAGN will be ab- 
sent, as its life-time is certainly short compared to the age 
of the Universe. But at present the DRAGN, which is sev- 
eral times larger than the cooling radius, is almost certainly 
depositing more energy into this region of the ICM than is 
being lost radiatively. 

Our best combined (PSPC and HRI) fit to the surface 
brightness profile, with a model consisting of a PSF, a /3 
model and a background, gives a /3 parameter of 0.74 ± 
0.03, typical for clusters of galaxies. The radio lobes are 
largely positioned beyond the X-ray core radius, allowing for 
projection, so they are expanding essentially into a power- 
law atmosphere with density falling as r"^'' ~ r~^'^^, quite 



close to th e r ^ profil e needed to give the lobes a self-similar 
structure iFalljll99ir) . 

The thermal pressure at the deprojected distance of the 
radio lobes is an order of magnitude larger than the mini- 
mum pressure of the lobes. Thus the minimum energy in the 
lobes is a severe underestimate of the actual energy content. 
This seems to be typical for DRAGNs, of both FR types, 
although only in a few cases can it be established as clearly 
as in Her A. 

Two features of the X-ray emission may result from the 
interaction of the cluster gas with the radio lobes. The re- 
gion around the core of the cluster is clearly elongated along 
the radio axis, as revealed by our on- and off-lobe profiles. 
In addition there are a number of apparently discrete X- 
ray enhancements projected on or around the radio lobes; 
however some of these may be artefacts or due to poor pho- 
ton statistics in our ROSAT data. Unlike some other well- 
studied DRAGNs, Her A as yet does not show X-ray holes 
coincident with the radio lobes. But we have shown that the 
such holes could be too shallow to detect in our data, due 
to projection effects: the main lobes being displaced along 
the line of sight away from the X-ray bright region around 
the cluster core. Deep observations with Chandra should re- 
veal some sign of these holes, if present, and also clarify the 
nature of the 'enhancements'. 
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